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Abstract — This paper develops upper bounds on the end-to-end 
transmission capacity of multi-hop wireless networks. Potential 
source-destination paths are dynamically selected from a pool of 
randomly located relays, from which a closed-form lower bound 
on the outage probability is derived in terms of the expected 
number of potential paths. This is in turn used to provide an 
upper bound on the number of successful transmissions that can 
occur per unit area, which is known as the transmission capacity. 
The upper bound results from assuming independence among the 
potential paths, and can be viewed as the maximum diversity 
case. A useful aspect of the upper bound is its simple form for 
an arbitrary-sized network, which allows insights into how the 
number of hops and other network parameters affect spatial 
throughput in the non-asymptotic regime. The outage probability 
analysis is then extended to account for retransmissions with a 
maximum number of allowed attempts. In contrast to prevailing 
wisdom, we show that predetermined routing (such as nearest- 
neighbor) is suboptimal, since more hops are not useful once 
the network is interference-limited. Our results also make clear 
that randomness in the location of relay sets and dynamically 
varying channel states is helpful in obtaining higher aggregate 
throughput, and that dynamic route selection should be used to 
exploit path diversity. 

Index Terms — multihop routing, transmission capacity, SINR, 
outage probability, stochastic geometry 

I. Introduction 

In a distributed wireless network with random node lo- 
cations, determining the precise network capacity is a long- 
standing open problem that includes many other simpler open 
problems as special cases IT]. Therefore, suboptimal analytical 
approaches that provide insights into the achievable throughput 
and inform improved protocol design are well-motivated, even 
if they fall short of strict upper bounds. Multihop routing 
is generally considered necessary in large wireless networks, 
both to ensure connectivity and to improve throughput, but it 
is typically not well-optimized, nor are its performance limits 
in general models known. In this paper we explore optimal 
multihop strategies by considering dynamic path selection. 
Pre-determined routing strategies such as nearest-neighbor 
routing, although they may perform fairly well on average, 
are generally not optimal for a given network state (which 
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includes node positions and all the channels among them). In 
fact, a randomly deployed set of potential relays with time- 
varying fading channels may allow a large gain by providing 
more potential routes. In this paper, we are interested in 
how the inherent randomness in the network can be better 
harvested to improve the end-to-end success probability and 
hence throughput over more static approaches. 

We assume that node locations are a realization of homoge- 
neous Poisson process in order to investigate the throughput 
statistically. This model provides a tractable way to charac- 
terize how the end-to-end success probability and throughput 
varies over different multihop routing strategies. We aim to 
study how multi-hop routing with the assistance of a pool of 
randomly deployed relays impacts the throughput scaling in a 
non-asymptotic regime, i.e. for networks of finite population 
and area, with the goal of finding how much the average 
throughput can be increased under quality of service (QoS) 
constraints. Considering uncoordinated routing selection, we 
aim at determining the fundamental limits for a general class of 
routing strategies instead of predetermined selection. It can be 
expected that the diversity gain resulting from the randomness 
and dynamic channels potentially provides significant through- 
put improvement. 

A. Related Work and Motivation 

The best-known metric for studying end-to-end network 
capacity is the transport capacity |2|-|j4]. This framework 
pioneered many notable studies on the limiting scaling be- 
havior of ad hoc networks with the number of nodes n 
by showing that the maximum transport capacity scales as 
Q {y/n) in arbitrary networks [2|. The feasibility of this 
throughput scaling has also been shown in random networks by 
relaying all information via crossing paths constructed through 
the network |5|. Several other researchers have extended this 
framework to more general operating regimes, e.g. ||6|, ||7). 
Their findings have shown that nearest neighbor multihop 
routing is order-optimal in the power-limited regime, while 
hopping across clusters with distributed multiple-input and 
multiple-output (MIMO) communication can achieve order- 
optimal throughput in bandwidth-limited and power-inefficient 
regimes. However, most of these results are shown and proven 
for asymptotically large networks, which may not accurately 
describe non-asymptotic conditions. Moreover, scaling laws 
do not provide much information on how other network 
parameters imposed by a specific transmission strategy affect 
the throughput. 
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If node locations are modeled as a homogeneous Poisson 
point process (HPPP), a number of results can be applied from 
stochastic geometry, e.g. |^|, |9|, in particular to compute out- 
age probability relative to a signal-to-interference-plus-noise 
ratio (SINR) threshold. These expressions can be inverted to 
give the maximum transmit intensity at a specified outage 
probability, which yields the transmission capacity of the 
network fTO |. This framework provides the maximum number 
of successful transmission the network can support while 
simultaneously meeting a network-wide QoS requirement. 
This framework allows closed-form expressions of achievable 
throughput to be derived in non-asymptotic regimes, which are 
useful in examining how various communication technique, 
channel models, and design parameters affect the aggregate 
throughput, e.g. |n|-|18|; see [19] for a summary. While 
the transmission capacity can often be expressed in closed- 
form without resorting to asymptotics, it is a single-hop or 
"snapshot" metric. Recent work |20|, [21] began to investi- 
gate the throughput scaling with two-hop opportunistic relay 
selection under different channel gain distribution and relay 
deployment. However, more general multi-hop capacity has 
not proven tractable. 

If several other strong assumptions are made, e.g. that all 
relays are placed equi-distant on a straight line and all outages 
are independent, then closed-form multi-hop transmission ca- 
pacity can be derived |22|. Stamatiou et. al. |23] also investi- 
gated multihop routing in a Poisson spatial model, whose focus 
is to characterize the end-to-end delay and stability, again 
based on predetermined routes. 'Other recent works analyzing 
the throughput of multihop networks using stochastic geo- 
metric tools include p4] , which extended p2| to non-slotted 
ALOHA and |251, which also adopted a similar framework to 
||25 1 to study the throughput-delay-reliability tradeoff with an 
ARQ protocol, and didn't require all hops to be equidistant. 
However, all of these used predetermined routing selection. In 
fact, the outage of a predetermined route does not preclude 
the possibility of successful communication over other routes. 
Separately, multihop capacity has also been studied in a line 
network without explicitly considering additional interference 
||26), p7| . This approach is helpful in comparing the impact 
of additional hops in bandwidth and power-limited networks, 
but fails to account for the interference inherent in a large 
wireless network. 

In addition, the above-mentioned diversity gain from dy- 
namic relay selection has been utilized for opportunistic rout- 
ing | [28) | p9] so any node that overhears packets can partic- 
ipate in forwarding. Reference f29l appeared to be the first 
investigation of the capacity improvement from opportunistic 
routing compared with predetermined routing in a Poisson 
field. However, the performance gain shown in f2S\ is based 
on simulation without an exact mathematical derivation. Dif- 
ferent random hop selection strategies have also been studied 
and compared | ,30j |31J without giving tractable throughput 
bounds. Hence, characterizing the available diversity gain is 
worth investigating. In this paper, we will explicitly show that 
since a pool of randomly located relays with varying channels 



provides more potential routes, more randomness is preferable. 

B. Contributions and Organization 

Instead of predetermined routing, dynamic route selection 
from random relay sets under varying channel states are inves- 
tigated in this work. The main contributions are summarized 
as follows. 

(1) We provide a lower bound on the end-to-end outage 
probability (Theorem [TJ, which can be expressed as an ex- 
ponential function with respect to the expected number of 
potential paths. This result implies that higher throughput can 
be achieved when the correlation among the states of different 
hops is low and hence randomness and opportunism is high. 

(2) We further derive in closed-form the expected number 
of S — V routes that can complete forwarding both for single 
transmission in each hop and for two different retransmission 
strategies that are subject to constraints on the number of 
allowed attempts, given in Lemma [T] The basic idea is to 
map all relay combinations to a higher dimensional space and 
focus on the level set with respect to the success probability 
function. 

(3) A closed-form upper bound on transmission capacity as 
a function of outage constraint e and the number of relays 
m for a general class of multi-hop routing strategies is given 
in Corollary [2j which follows from Lemma [T] and is the main 
technical result in the paper. These closed-form results assume 
a general exponential form of success probability, which 
includes most commonly used channel models as special cases, 
including path loss, path loss with Rayleigh fading, and path 
loss with Nakagami fading. 

The above results show that in networks with uncoordinated 
routing, an "ideal" diversity gain arising from independent 
statistics of different paths allows the throughput to exhibit 
near linear scaling in the number of relays m as long as 
the density of relay nodes exceeds the threshold imposed by 
the outage constraint e. This diversity gain requires strong 
"incoherence" among different paths, which would presumably 
degrade for large m since longer routes are more likely 
to be correlated or share common links. Unlike the single- 
hop scenario where network throughput must decrease about 
linearly as the output constraint e is tightened, the multihop 
capacity bound is less sensitive to e especially for large m. 

(4) Finally, we briefly show that all predetermined routing 
strategies with no central coordination and without further 
information like channel state information (CSI) may fail to 
outperform single hop transmission in an interference-limited 
network because of the large increase in interference. Hence, 
exploiting randomness is important for multihopping to be 
viable in networks of finite size. 

The rest of the paper is organized as follows. In Section 
[m we first define the end-to-end metric that quantifies the 
network-wide throughput, and then state the key assumptions 
for the analysis, as well as the channel models and their 
associated general form of per-hop success probability. We 
then develop and derive lower bounds for end-to-end outage 
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probability for general scenarios in Section |III-A| Specifically, 
this provides a closed-form lower bound if the channel model 
allows the per-hop success probability to be expressed in 



exponential form, which is developed in Section III-B This in 
turn results in an upper bound for the multi-hop transmission 



capacity in Section III-C Implications and interpretations of 
the results are provided in Section |IV] 

II. Models and Preliminaries 

A. Models and Assumptions 

We assume that the locations of all sources are a realization 
of an HPPP Et of intensity Aj, and a set of relays are also 
randomly deployed in the plane with homogeneous Poisson 
distribution independent of Sf. We consider a fixed-portion 
model, i.e., the relay set is of spatial density ^^-^Xt, where 
7 G (0, 1) is assumed to be a fixed constant. In other 
words, if the locations of all wireless nodes are assumed 
to be an HPPP S with intensity A, then the set of active 
transmitters is of intensity At = 7A. The destination node is 
assumed to be a distance R away from its associated source 
node in a random direction, and is not a part of the HPPP. 
Suppose transmission rate 6 ~ 5 log (1 + f3) is required for 
successful transmission, where /3 is therefore the required 
SINR. Denoting e as the target outage probability relative to /?, 
the transmission capacity [lOJ in an uncoordinated single-hop 
setting is defined as 



T(e) = (1 - e) max At 

P(SINR</3)<e 



(1) 



which is the maximum expected throughput per unit area. 
Since b is simply a constant function of (3, we ignore it for 
simplicity. 

Now suppose each session uses k transmissions with the 
assistance of the relay set. These k attempts can be performed 
in an arbitrary k orthogonal slots, i.e., the unit time slot can 
be divided into k equal subslots and the source and relays 
take turns transmitting in these subslots: only one transmitter 
per route is active at a time. The contention density is still of 
density At, but each packet is transmitted k times. Therefore, 
the multihop transmission capacity metric should be modified 
to be 



max , 

P(SIR<^)<e k 



(2) 



since each hop requires a time slot, so the overall throughput 
must be normalized by k. It should be noted that although one 
can "pipeline" by simultaneously transmitting different packets 
on different hops, this does not change the transmission 
capacity since the contention density simply becomes k times 
larger Similar analysis can be applied in quantifying the 
transmission capacity with this intra-route spatial reuse but 
leads to the same result. When no retransmissions are allowed, 
we have k ^ m + 1 with m relays; if we consider M total 
attempts (including retransmissions) for any single session, 
then k ^ M. 

Slivnyak's theorem p2| states that an entire homogeneous 
network can be characterized by a typical single transmission. 
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Table I 

Summary of Notation and Parameters 
path loss exponent 

SINR requirement for successful reception per hop 
outage probability constraint 
S — T) transmit distance 
contention density of all potential transmitters, 
intensity of H 

contention density of active transmitters in any 
subslot. intensity of Ht 
the portion of active transmitters. At = A7 
Euclidean distance from node i to node j 
general parameters in exponential form single-hop 
success probability, 
90 (rii,M) = Gcxp (^-MKr'f^^ 

end-to-end outage probability employing m relays 
without retransmissions 
Zm the corresponding location vector in 

90 (Xij,Xt) single hop success probability 

9m (-Zmi At) the end-to-end (m + 1 hops) success probability 

with the relay set Z,„ 
K, G general parameters in exponential form single-hop 

success probability, 
90 {rij,Xt) = Gcxp (^-XtKr^.^ 
f2m (B) Lebesgue measure of B 

dm{Zm) sum of squared distance of m + 1 hops with the 

relay set Zm 

Dm, distance constraint, dm (Zm) < Dm 

Tm (e) m + 1 hop transmission capacity 

A A := X-^K 

K K := Gtt{1 — 'y)/-fK 



Conditioning on a typical pair, the spatial point process is 
still homogeneous with the same statistics. Suppose that all 
transmitters employ equal amounts of power, and the network 
is interference-limited, i.e. noise power is negligible compared 
to interference power. Relays can be selected from all nodes 
in the feasible region. In this paper, we consider the effects 
of both path loss and fading. For point-to-point transmission 
from node i to node j at a distance Vij, the requirement for 
successful reception in this hop is expressed in terms of signal- 
to-interference ratio (SIR) constraint as 

2 



SIR^ 



> A 



(3) 



Efc^^ill^feill '-kj 

where a denotes the path loss exponent, and hij the fading 
factor experienced by the path from i to j. Distinct links 
are assumed to experience i.i.d. fading, which is typically 
reasonable. 

B. Per-Hop Success Probability 

A Poisson node distribution often results in an exact or ap- 
proximate exponential form for per-hop successful probability. 
That is, given that the packet is transmitted from node i to next 
hop receiver j over distance r^j and contention density At, the 
probability that the received SIR stays above the target /3 can 
be expressed as 



P(SINR,j > /?) = 9o{n,,Xt) = Gexp (-Ati^r^) 



(4) 



where G and K depend on the specific channel models and are 
independent of r^j and At. This holds for several commonly 
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used channel models, including Rayleigh fading with path loss, 
Nakagami fading with path loss, and path loss without fading, 
as we briefly show in this subsection before using the general 
form in the remainder of this paper 

1) Rayleigh fading: Baccelli et. al. fSl showed under 
Rayleigh fading that: 



go{rij,\t) = exp 



(5) 



where C{a) = 27rr(f)r(l - with V{z) = 

Jp°° exp(— t)dt being the Gamma function. Hence, the 
coefficients under Rayleigh fading can be given as 



/^RF = /3°C(a); Grf = 1. 



(6) 



2) Nakagami fading: Nakagami fading is a more general 
fading distribution, whose power distribution can be expressed 
in terms of fading parameter mo > 0.5 as 



fziz) 



r(TOo) 



■ exp(— TOoz) 



(7) 



Recent work fTSl has suggested a way to study the outage 
probability by looking at the Laplacian transform. In both low- 
outage and high-outage regimes, the success probability under 
Nakagami fading can be expressed as an exponential function 
with the following coefficients 

low outage regime: Kj^^ ~ rim„/3° , Gnf = 1, 
high outage regime: i^NF = f^mo/?", 



fe=i ;=i 



fc! 



Both the derivation of these coefficients and the definition of 
and T^ j can be found in Appendix [A| These two regimes 
are typical in practical systems. 

3) Path Loss Model {non-fading): The exact closed-form 
formula of the success probability with only path loss is 
unknown. One approach is letting toq ^ co for Nakagami 
fading, which converges to a path-loss-only model. However, 
lower and upper bounds follow an exponential form. For 
example, partitioning the set of interferers into dominating 
and non-dominating nodes, an upper bound can be obtained 
as g^{rij,\t) = exp(— At7r/3a r^^). Reference flO^ has also 
shown an upper bound on the transmission capacity that is 
times the lower bound for small e, and has illustrated by 
simulation the tightness of these bounds. This suggests that 
there exists some constant ifpL such that 



a^irij, Xt) = exp (-AtifpLry), 



(8) 



where 



2 cy. 2 

^/3=<ifpL< GpL = l. (9) 

a — I 



III. Main Results 

In this section, we first address the end-to-end outage prob- 
ability for general channel models, which builds a connection 
between outage probability and expected number of potential 
paths. Next, if channel models allows the per-hop success 
probability to be expressed in an exponential form, closed- 
form lower bounds can be derived. This in turn will provide a 
closed-form upper bound for multihop transmission capacity. 

A. Outage Probability Analysis for General Per-hop Success 
Probability 

Suppose that m relays are employed by a typical source- 
destination iS—T)) pair. Since all S—T) pairs are stochastically 
equivalent, we can investigate the performance by looking at 
a typical S — T> pair. We will build a connection between 
the outage probability and the expected number of relay sets 
that can connect the source and destination. Suppose that there 
is a transmission pair with source and destination located at 
(-i?/2,0) and (i?/2,0), respectively. With the < i < 

m) relay located at {xi,yi), let Zm = (xi, . . . , a;„, y™) 
denote the location vector of this specific relay set. From 
Slivnyak's theorem, conditional on a typical transmission pair 
or finite number of nodes, the remaining point process is still 
homogeneous Poisson process with the same spatial density 
(we ignore a finite number of singular points here). Therefore, 
all relay combinations form a homogeneous point process in 
a 2m-dimensional space 7?.^™, as illustrated in Fig. [l] The 
effective spatial density is A = (1 — 7) A, which characterizes 
the density of the pool of nodes that have not been designated 
for specific transmissions. Assume that each relay combination 
Z„i can successfully assist in communication between the 
S — T) pair with probability g^iZm, At). If we call a relay 
set that can successfully complete forwarding in a given 
realiazation of the spatial process a potential relay set, then 
the expected number of potential relay sets in a hypercube B, 
denoted by Nb, can be expressed as 

lim E(Nb) 
^ lim E ^ I (Z„, is a potential relay set) 



(10) 



where V2m{B) denotes the Lebesgue measure of B and I(-) 
denotes indicator function. Let the random variable iV^ be 
the number of relay sets that can complete forwarding using 
m + 1 hops. Then E(iVm) is the expected number of different 
routes that can successfully forward the packets for an 5 — I? 
pair. A larger E(iV„i) naturally leads to lower outage, where 
we formalize in the following theorem. 

Theorem 1. Assume that all end-to-end transmissions are 
achieved via m + 1 hops with m relays. The end-to-end outage 
probability for any S — T) pair p^^^ can be lower bounded as 

pir^^ > exp(-E(7V„)). (11) 
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Figure 1. In the left plot, the S—T) pair use 3 relays to assist the transmission, 
which can be matched to a point in high-dimensional space as plotted in the 
right plot. In fact, a precise plot requires drawing on 6 dimensional space with 
each relay accounting for 2 dimensions, but the right plot may help explain 
the mapping intuitively. 



Proof: The key idea is to view the outage event as the 
intersection of a set of decreasing events. The basic properties 
of decreasing events suggests a lower bound by treating all 
these events as mutually independent. See Appendix [B] for 
complete proof. ■ 

The lower bound can only be approached when the pool 
of potential relay sets form a Poisson point process in a 
corresponding 2m-dimensional space (detailed in Appendix 
[B] i, i.e. all potential relay combinations independent from 
each other. In practice, route selection for different source- 
destination pairs are not independent, so the derived lower 
bound is not obtained by realistic routing strategies. This 
result, however, indicates that low correlation among differ- 
ent routes can reduce the outage probability in essence by 
enhancing diversity. 

We conjecture that this bound is tight and reasonable for 
small m (e.g. the bound is exact for single relay case) but 
may be loose for large to. This is because for a fixed pool of 
relays, the correlation among different routes increases when 
the number of relays m grows, i.e. for large to, many routes 
are likely to share at least one link. Also, the bound may 
become loose for low outage regime (e.g. e < 10""^), because 
the outage bound is an exponential function of E (iV,„), where 
even constant factor difference may result in an exponential 
gap. The advantage of this bound is to allow us to see how 
the maximum allowable throughput scales with fixed normal 
outage constraint in the non-asymptotic regime, which will be 



shown later In fact, some constant correction factor can also 
be applied based on some simulation results without changing 
the scaling of maximum contention density versus outage 
constraint. The expected number of different routes E [Nm] 
plays an important role, which will be calculated exactly in 
the following subsection. 

We caution, however, that the key assumption of indepen- 
dent path selection cannot be achieved in practice. Potential 
routes are virtually coupled and never independent from each 
other For instance, failure of finding a potential route for 
one session typically implies a lower success probability for 
another session. But when uncoordinated routing selection 
is employed for all S — T) pairs and when the number of 
relays is reasonably small, the correlation is mitigated. And 
this simplication assumption allows a reasonable closed-form 
bound to be derived. 

B. Outage Probability Analysis for Exponential-Form Per-hop 
Success Probability 

Now we begin to concentrate on the success probability 
of exponential forms. When no retransmissions are adopted, 
if a specific route is selected for packet delivering over to 
relays with hop distances ri, r2, ■ ■ ■ , ^m+i, respectively, 
the probability for successful reception can be found as the 
product of each hop's success probability: 

, Tm+l, At) 



exp 



m+l ^ 



(12) 
(13) 



m+l 

i=l 

Here, we assume independence among the success probability 
of each hop, which requires sufficient diversity in the interferer 
locations over time. Although this assumption is not valid 
in practice since transmissions across subsequent hops are 
correlated, it allows us to retain tractability and has been 
shown to be a reasonable approximation in uncoordinated 



networks (as illustrated in |22|). The temporal correlation may 
also be mitigated through diversity techniques like frequency 
hopping. 

Conditional on a typical transmission pair with source and 
destination located at (— i?/2,0) and (i?/2,0), the spatial 
point process is still a HPPP with the same statistics. In the 
TO relay case with the < i < m) relay located at 

(xi.yi), let Z„ = {xi,yi,. . . ,Xm,y7n) denote the locations 
of the specific relay set, then we can define the corresponding 
distance statistics as 



A R 

d,n{Z,n) ={xi + — f + 



m — 1 
Y{X^ 
1=1 
-1 



\ 2 / ^ \ 2 

Xi+i) + [Xrn - -^j 



1=1 



(14) 



This is the sum of squares of hop distances. Hence, the routing 
success probability for a specific set of relays with location 
Zm can be explicitly expressed as 

gm{Zm,\t) = G""+' exp(-Atifd,„(Z„)). (15) 
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In fact, an arbitrary set of relays will have positive probability 
for successful forwarding. However, for those relay sets with 
large d„i{Z„i), the communication process becomes extremely 
fragile and difficult to maintain due to the low reception prob- 
ability and large distance. Practical protocols usually attempt 
to search potential routes inside a locally finite area instead of 
from the infinite space since the longer routes are very unlikely 
to be an efficient one. In order to leave the analysis general, 
we impose a constraint dm{Zm) < Dm for the m relay case, 
where Dm — > oo reverts to the unconstrained distance case. 
We will later show that a reasonably small constraint Dm is 
sufficient to achieve an aggregate rate arbitrarily close to the 
capacity upper bound. 

Moreover, since only one transmitter is active at a time 
along the entire multihop route, in each subslot, each node 
is used as a relay by other source-destination pairs with 
probability 7. Therefore, the pool of relays in each hop can be 
treated as the original point process S with each point being 
deleted with probability 7. Hence, the location of all relay 
sets in M"^™ can be viewed as a realization of a point process 
with effective spatial density A™(1 — 7)™. This leads to the 
following lemma. 

Lemma 1. Define n = G7r(l — "f)/jK and A = X^fK. If all 
end-to-end transmissions are achieved via m + 1 hops with m 
relays, with a constraint dm{Zm) ^ Dm (Zm G ^^'"), the 
expected number of potential relay sets can be computed as 



m + 1 

m— 1 



m + 1 



exp(-AAn) 



771+1 



(16) 



Proof: The key point in the proof is that the isosurface 
of dm{Zm) forms a high-dimensional elliptical surface, which 
provides a tractable closed-form solution. See Appendix |C] ■ 
This result indicates that a larger m typically provides more 
diversity, because it provides more possible combinations of 
different relays, and the dynamically changing channel states 
provide more opportunities for us to find a potential route. 
A larger feasible range for route selection Dm also increases 
the expectation, but since the effect of Dm mainly exhibits 
as an exponentially vanishing term, it can be expected that 
a fairly small range is enough to approach the limits. More- 
over, this analytic framework can be extended to account for 
retransmissions in the following two scenarios. First, the best- 
effort retransmission protocol requires that each hop adopts k 
retransmissions regardless of the results of each transmissions. 
The following lemma provides more general results for best- 
effort protocols by allowing each hop to adopt a different 
number of retransmissions. Second, instead of specifying 
retransmissions for each hop, we bound the maximum number 
of total allowed attempts to M. Define 1 := (1, • • • , 1)"^. The 
following lemma provides closed-form results for these two 
scenarios. 



Lemma 2. Assume that all end-to-end transmissions are 
achieved via to + 1 hops with m relays. 

(1) In the best effort retransmission setting, if the < 
i < m + 1) hop is retransmitted ki times, then the expected 
number of potential relay sets can be given as 



Ek(A^™) 



E 



(-i)™+V'"(i-7r(-G) 



2^i = l 



n: l^n^k 

m+1 



n 



k, \ "^^"H'^rf^ 



(17) 



(2) If the S — T) transmission allows M transmissions in 
total without specifying the number of retransmissions for each 
hop, then the expected number of potential relay sets can be 
given as 



m+l 



31 - 1 



k^ l<M j^ l<A/ ' = 1 

k^l j^k 

(18) 



where Ek {Nm) is given in \17\ . 

Proof: See Appendix [D] ■ 
Similarly, this lemma is derived by mapping all potential 
relay sets onto a 2TO-dimensional space and investigating the 
isosurfaces in that space. The above results on the expected 
number of potential relay sets immediately yield the following 
corollary. 

Corollary 1. Assume that all end-to-end transmissions are 
achieved via to + 1 hops with m relays. 

(1) If only a single transmission is allowed for in 
each hop, then the outage probability under a constraint 
dm{Zm)l^Dm (Zm^^"^"^) Can be computed as 



Gk^ 



exp 



TO+1[_ V TO+1 

m— 1 



Ai?2 



exp {-KDm) E ^ ( A [D, 




(19) 



(2) If each hop adopts kiil < i < m -\- 1) retransmissions, 
the outage probability can be computed as 



> exp - ^ 



n: l-<n-<k 



(-1)'" + ^ TT" (1 - 7)" i-G)^T=V ■ 



i=l 



AR-' 



(20) 
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(3) If the transmission adopts M retransmission in total, 
the outage probability qm can be given as 



( 



qM > exp 



E n 



ki-l 



Ek (iV™) , (21) 



where Ek (Nm) is given in (17 \. 



This corollary provides closed-form lower bounds on the 
end-to-end outage probability. For sufficiently large D„i in 
the non-retransmission case, the lower bound reduces to: 



(m) ^ 

Pout > exp 



exp 



j"^K"'{m + 1) 



TO + 1 



(22) 



which gives a clear characterization for low-coherence routing 
selections. As expected, multi-hop routing with the assistance 
of randomly deployed relays improves the success probability 
by providing large potential diversity, with the randomness in 
both relay locations and channel states proving helpful. 

We note that unlike the single hop scenario |10|, our bound 
for outage probability without retransmissions is not globally 
monotonically increasing with A if I?,,„^oo. For sufficiently 
large but not infinite D^, the outage probability can be 
approximated through a first-order Taylor expansion in the low 
density regime: 



> exp 
« 1 - A 



m + 



TO 



1 



exp 



Dm 



AR^ 

TO + 
2 



exp (- A An) 



R 



1 



(23) 



which indicates large outage probability in the low density 
regime, arising from the difficulty in guaranteeing a relay 
within range dm{Zm) < Dm in a sparse network. The detailed 
monotonicity can be more closely examined by studying the 
function /(A) = exp(— aA) — exp(— 6A) {b > a > 0), whose 
derivative can be computed as: 



/'(A) =exp(-6A){6-aexp[(6-a)A]}. 



(24) 



The maximum value of /(A) occurs at Aq = 53^ In ^, and 
/(A) is monotonically increasing at (0, Aq] and decreasing at 

we 



(Ao, 00). Using this property, and defining A — (^^+1)0 
can see that 



(m) / \ \ ^ 



, A 

Al-A 



where the minimizing A is: 



lK[Dm- 



m+l 



In 



(to + 



(25) 



(26) 



Hence, pilJ^''(A) is monotone in both [0, Ao] and (Ao,oo). 
Taking the inverse over (Ao,oo) will yield the bounds on 
maximum contention density. 



C. Transmission Capacity Upper Bound 

When Dfn — > 00, Aq goes to 0. Therefore, Pom\^) is 
monotonically increasing in (0,oo). Therefore, we can get the 
following transmission capacity bound by taking the inverse 
of this outage probability function. 

Corollary 2. (1) If each hop adopts a single transmission, the 
transmission capacity can be bounded as 



TO In ^^^^^ + InG - ln(m + 1) - In In i 



KR^ 



(l~e) 
(27) 



where e > exp (-^^^ 

(2) If best effort retransmissions is adopted with each hop 
utilizing ki retransmissions, the transmission capacity can be 
bounded as 



,m+l > 



an^ 



In {Gm) 



-Inlni 



(28) 



where e > exp 



Proof: See Appendix |E] ■ 
In the case where each hop use a single transmission, when 
Dm 00 but is reasonably large, the outage probability can 
be approximated using L' Hospital's rule 



e > exp 



1 



exp 



KR^ 

TO + 1 



exp 



TO + 1 



By simple manipulation, the upper bound T,"''(e, L)„i) on 
transmission capacity with Dm constraint becomes 



T,"„\e, Dm) + e{exp{-T;;''(e, Dm)K{m + l)An}} 



(29) 



which means the gap between the general bound and the bound 
with distance constraints will decay exponentially fast with 

Dm- 

It should be noted that if we only impose a constraint on the 
maximum number of allowable attempts A/, it is difficult to 
get a closed-form capacity bound. But since the outage bound 
is monotonically increasing with A, it would allow a numerical 
solution. 



Figure 2. The end-to-end outage probability pout computed numerically and 
its lower bound versus the contention density A of all potential transmitters. On 
each data point on the simulation curves, an error bar indicating the confidence 
interval (which is chosen to be the emperical standard deviation) is plotted. 
Note that the lower bound is fairly tight for transmissions with 2 hops and 3 
hops, and becomes looser for transmissions through 4 hops. 



IV. Numerical Analysis and Discussion 

In this section we study the implications of the theoretical 
results through simple numerical analysis and simulation. 
The presented plots presume path-loss attenuation, Rayleigh 
fading, and no noise. The SIR threshold f3 is set to 1 while 
other parameters are varied. The S — T) distance is primarily 
R — A although i? = 6 is also used. In a interference-limited 
environment the exact values of R and A are not particularly 
important since the outage probability is constant for fixed 
\R^. In the simulation, a spatial Poisson point process is 
generated. For each spatial density A, we pick S — V pairs 
uniformly at random that has an average spatial density A7, 
and let each S — T) pair perform uncoordinated path selection. 
If there is any hop conflicting with routes selected by other 
S — D pair, then the transmission fails. The set of connected 
paths from 5 to 2? is determined through SIR measurements of 
each link, and constrained to a maximum end-to-end distance 
of Dm = 60^ which qualitatively approximates the extreme 
case of D„j. 00 for R < 6. 



Tightness of outage lower bound. The lower bound (22 1 is 
plotted against simulated outage probability in Fig. [2] The sim- 
ulated outage probability takes into account the dependence 
among consecutive transmissions and parallel path selection. 
For each data point, an error bar is plotted to indicate the 
confidence interval of the simulation results. Here, the width 
of the confidence interval is chosen to be twice the empirical 
standard deviation. The bound is observed to indeed be a lower 
bound and to be quite tight, albeit slightly looser for increasing 
numbers of hops. 

The Number of Relays m. Since 9(ln(m + 1)) is neg- 
ligible compared to 8(m), the transmission capacity bound 



Figure 3. The maximum allowable effective contention density A/(m + 1) 
computed numerically versus the number of hops when i? = 4, o = 3, /3 = 1 
for a typical S ~ T> pair It can be seen that the maximum effective density 
scales nearly linearly in the number of hops. Note that increasing the distance 
constraint D only provides fairly small throughput gain. 




■ ^ - 3 hops 
— H — 4 hops 
-H — ■ 5 hops 
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Figure 4. The maximum allowable contention density A computed numeri- 
cally versus the intensity ratio of relays to source nodes 
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when e = 0.05 



(27 1 exhibits near linear scaling behavior with respect to the 



for a typical S — T) pair. It can be observed that the maximum density scales 
logarithmically in 



number of relays m. This gain arises from the increasing 
route diversity as m grows, since more hops allows more 
potentially successful routes. This gain does not depend on 
the noise level and is not achieved by pre-determined routing 
approaches, which primarily are useful for overcoming per-hop 
range limitations (i.e. noise). We caution that this upper bound 
is likely to be increasingly optimistic for large m, since longer 
potential routes will presumably result in higher correlation 
between candidate paths. Fig.[3]shows the maximum allowable 
contention density versus the number of hops for different 
outage constraints. As expected, the effective contention den- 
sity scales nearly linearly for small m as expected, and then 
diminishes rapidly for large m. In practice, a modest number 
of hops would be taken since longer routes experience larger 
delay and more protocol overhead. The proper choice of m 
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under realistic correlation and protocol overhead models is an 
interesting topic for future research. 

Outage Probability Constraint e. The transmission ca- 
pacity bound is not sensitive to the outage constraint e in 
the low outage regime, because the double logarithm as in 
In In - largely reduces its sensitivity. For instance, when the 
target e is decreased from 10^^ to 10^"*, the throughput 
only experiences a small constant loss. This is quite different 
than single-hop transmission capacity, which exhibits linear 
scaling with e in the low outage regime and so going from 
e = 10~^ to e = 10^^ would in fact decrease the transmission 
capacity by two orders of magnitude [lOj. Hence, multihop 
transmission capacity is apparently much more robust to severe 
QoS constraints compared to single hop. 

Availability of Relays. Recall that nodes in the network 
are divided into a fraction 7 that may transmit and 1 — 7 that 
are available as relays. Corollary [2] implies that increasing the 
pool of relay nodes will logarithmically increase throughput, 
so the diversity gains diminish rapidly once a large enough 
pool to guarantee multi-hop route selection exists. Note that 
we primarily consider fixed-portion relay models here, which 
means the density of the pool of relays grows along with 
the density of source nodes. Simulations in Fig. |4] show the 
maximum contention density versus the intensity ratio ^^-y- 
of relays to source nodes, with e = 0.05. The results can 
be modified to study fixed-density relay models (where the 
density of relay nodes is a fixed constant Aj.) by substituting 
with Xr/Xt, which we do not present here. 

Sum-Squared-Distance Constraint D,„. The gap between 
the distance-constrained maximum density Tm(e, D„i) and the 
transmission capacity Tm(e) is subject to exponential decay 
with respect to (m + 1)1?™ as predicted in (29 1. Hence, 



searching for multihop routes in a local region should be 
sufficient. Fig. [3] illustrates this when the S — T) distance 
is i? = 4. It can be observed that when D is reasonably 
large compared with i?^, increasing provides almost no 
throughput gain. Also, this gain shrinks rapidly as m increases. 



which can also be expected from (29 1. 

Limitations and Future Directions. The results of this 
paper are well-suited to both fading and non-fading channels, 
but care should be exercised in considering more diverse 
channel models like log-normal shadowing, which do not nec- 
essarily lead to an exponential outage probability expression 
|[TT|. In addition, the models in this paper assume mutual 
independence among different links, which would not hold in 
general, particularly for routing and scheduling strategies that 
require cooperation. Furthermore, the theoretical gap between 
our upper bound and the true capacity is unknown, and how 
it grows with the number of hops is of interest. 

The multihop transmission capacity bound shows that dy- 
namic routing selection is of significant importance when 
there is sufficient randomness in the network as far as path 
gains and interference. In fact, predetermined routing (like 
nearest-neighbor) is unlikely to approach the throughput bound 
in interference-limited networks. A simple argument shows 
this. Considering a typical source-destination pair, the outage 



probability can be bounded as 

1 -pim^ = Gexp(-A7ifd„(ZM)) 

i?2 



< Gexp ~X-fK 



1 



(30) 



The equality can be achieved if and only if the m relays are 
equally spaced along the line segment between source and 
destination. In fact, from the properties of Poisson random 
process, this is almost surely unlikely to occur, resulting in a 
strict inequality. Setting A7(l — e)/(m+ 1) to r^(e), we can 
immediately get an upper bound 



1-e, G 
in - 



(31) 



KR'^ ~" l-t 

which is exactly equal to the single hop case. This suggests 
that predetermined routing will not provide further throughput 
gain in interference-limited networks compared with single 
hop direct transmission. 

We note that the power-limited regime (i.e. including noise) 
is not considered in this paper. Although noise is unimportant 
in the high density regime, it can be quite important in the low 
density regime, which is often power-limited. Our framework 
is primarily based upon an exponential form of per-hop success 
probability, which does not hold in the low SNR case. From a 
capacity perspective, the high density case is of more interest, 
since in the low density (power-limited) regime, nodes have 
far fewer options as far as selecting relays, and spatial reuse 
in the network is not very important. In fact, multihopping is 
known to be particularly helpful in changing a power-limited 
network to an interference-limited one by increasing the SNR 
in each hop, consistent with |22) , p6) . 

The design of transmission strategies that exploits the path 
diversity gain are left for future work. We conjecture that 
hop-by-hop route selection — which is much more realistic 
in a distributed network than the complete route selection 
assumed here — will achieve a lower diversity order (and 
hence transmission capacity). 

Appendix A 

Single hop success probability under Nakagami 

FADING 

The single-hop success probability with Nakagami fading 
can be developed as 



go{rij,Xt) 



E 

fc=0 



(32) 



where Ci^ (s) is the Laplace transform of the general Poisson 
shot noise process, and £^^■'(5) denotes the /c* derivative of 
£/^(s). The closed-form formulas of them are given by 1 15| 
as 



>C/^(s) = exp < -Atfi,, 



(33) 
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exp <-Xt^r, 



\mo J 



E 

3 = 1 



mo) 



2-1 ] 



(34) 



where Tkj is a constant defined in | [T5) , and 

rno- 



o "'0-1 , 

iTT \ ^ / m 



fc=0 



, ,B{k+-,mn-k--) (35) 



region Ti, we have — 0"=! Consider the hypercube 
Ti as . . .,Xm,ym), (xi + Sxi,yi + 5yi, . . . ,x,„ + 

Sx,n,ym + Sym)] when (5xj— >-0 and Jy^— >-0. Define = 
{xi,yi, . . . ,Xm,ym)- Since this is a simple point process, 
we can approximate the void probability as follows if the 
Lebesgue measure V2m(J'i) is small or n is sufficiently large 

lim ¥{Ai) 



lim 1 - g,n{Zi, At) TT ( 1 - cxp{-XSxiSyi) 



with B{a,b) denoting Beta function. By manipulation, we 
have 



1-1 fe 



x-^ v-^ 1 



A;=l 1 = 1 



- k.l 



lim I- g,n{Zi,Xt)W~X8xi5yi 
lim exp ( -A"5m(Zi, Af)u2m(-7^i) 

oXi=^Q,5yi~^0 \ 



(37) 



(38) 



Generally speaking, this does not have an expected expo- 
nential form. But we can simplify the expression in certain 
cases. For small single-hop outage constraint e, we have 
At^^mo/^" therefore go{rij,Xt) can be approximated 



as 



5o(?'y , At) « exp |-Atf],„„/3°ry I 



Consider two realizations lo and uj' of this higher-dimensional 
point process, and denote u < lo' if uj' can be obtained from 
ijj by adding points. An event Ai is said to be increasing 
if for every uj < lo' , (w) < I^. [uj') with I^. denoting 
the indicator function of the event Ai. If Ai(l<i<n) are all 
increasing events, then the Harris-FKG inequality p3) yields 



(36) 



In contrast, for large single-hop outage regime, i.e. 
Xt^moP^ffj ^ 1' employing L'Hospital's rule yields 



P(^)-P(f|A)>niP(A:). 



(39) 



9o[rij,Xt) 

( niQ — l k 



Letting n go to infinity, we can get the lower bound of outage 
probability as follows 



I, k=i 1=1 ^ 
We summarize them as follows 



exp|-Atano/3°'^*jj- V{A)> \h-aT\¥{Ai) 



low outage regime: i^Np — ilmoP" , Gnf 
high outage regime: i^NP = f^mo/?", 



1=1 



exp -X'"g^{Zi,Xt) lim V^U2m(-7^j! 

\ n— foo 



i=l 



exp(-E(A^„)). 



(40) 



fc! 



'^k,l 



fe=l 1 = 1 

Since practical system typically require low outage probability, 
our analysis may still work to a certain extent. 

Appendix B 
Proof of Theorem[T] 

Let the high-dimensional feasible region T for relay sets 
be the allowable range to select relays from determined by 
different routing protocols and design parameters. Denote by 
A the event that there is no relay set within T that can 
successfully complete forwarding. Ignoring the edge effect, 
we attempt to approximately divide T into n disjoint hy- 
percubes ^i(l<i<n) each of equal volume. For sufficiently 
large n, this approximation is exact. Let Ai{l<i<n) be the 
event that there exists no potential relay set within that 
can complete forwarding. Since the outage event A occurs 
only when there is no potential relay set in any of the 



Appendix C 
Proof of Lemma[T] 

The isosurface of = a has the following coordinate 

geometry form 



(41) 



where 



2 m-l 



-'^sum = ( 2^1 + ) + E ~ + ' ^ ^ 



Y,um = yl + ^{y,+i-yif + yli 



If we treat Xi^yi (1 < i < m) as mutually orthogonal 
coordinates, then ( |4T] ) forms a quadratic surface in 2m- 
dimensional space. See Fig.|2]for an illustration when m = 1. 
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From the properties of quadratic forms, the x part and y part 
of ( pT) can be expressed as 



Y 

-L SI] 



(CX - Rx)'^Ax(CX - 
(CY)TAy(CY), 



(42) 



where C, C are orthogonal matrices, Ax, Ay are diagonal 
matrices, Rx is a rn-dimensional vector, and is a constant 
that will be determined in the sequel. Here, the orthogonal 
transformation of X (Y) by C (C) and translation transfor- 
mation by Rx only result in rotation, flipping or translation of 
the quadratic surface without changing the shape of it. Since 
the corresponding quadratic terms of Xjum and Ysum have 
equivalent coefficients, we have A„j = Ax = Ay. Denote 
the symmetric quadratic-form matrix corresponding to ysum as 
Am, then A,„ is the following tridiagonal matrix of dimension 
m: 



A 



/ 2 -1 

-1 2 -1 
0-12 



V 










2/ 



(43) 



In fact, A„i is the canonical form of A„j with its eigenvalues 
on the main diagonal. Through orthogonal transformation and 
translation, X^um, Ysum can be brought to the explicit form: 



Y 



E 

1=1 



i=l 



XiXi 



Km 



^ 2 



(44) 



where Xi^iji are the new orthogonal coordinates and is 
the zth eigenvalue of A^. By its definition, Xjum is positive 
definite, and the following minimum value can be obtained 
if and only if m relays are placed equidistant along the line 
segment between the source and destination: 



> 



{\^+Xi\ + \x2~Xi\ + --- + \x„,- 



> 



m 



Therefore, i,; 



1 

r?j.+l ■ 



Now, (|41|i can be brought to: 

i?2 



hiji^ = a- 



1 



(45) 



(46) 



From the positive definiteness of A„j, A, > holds for all i, 
i.e. the above equation forms the surface of a 2TO-dimensional 
ellipsoid. Fig. |5] illustrates the ellipsoid when to = 1, which 
reduces to a circle. The Lebesgue measure of the ellipsoid can 
be written as: 



m+1 ' 



-Y 



!nr=iA. 



to! det (A„ 



(47) 




Figure 5. In single relay scenario, the contour of end-to-end success 
probability gi = cxp (— AtX (r^ +^2)) exactly a circle. The plot 
is a realization of Poisson point process of A = 0.1 /unit area. The source S 
and destination D are a distance 30 apart with relay R on the dotted circle 



satisfying -(- = 500. 



We also need to determine det (A™), which can be computed 
by the Laplace expansion of the determinant 



det(A,„) = 2det(A„j_i) - dct(A„_2). 



(48) 



Solving this recursive form with the initial value det(Ai) — 2 
and det(A2) — 3 yields 



det(Am) = TO + 1 



Kn(a)^ ,":+^^ . (49) 



(m -I- 1)! 

Now, we can compute the outage probability. Integrating 
over different isosurfaces with g{Zjn,\^) = exp(— A7iira), 
and defining H = ^'^'^k ^ ' can compute the average 
number of potential relay sets as 



A'"(l-7)' 



''^^'"^''^-G"+iexp(-Aa)da 



G 



'1+1 

m+l \ ^11 



A"(l - 7)" 



TO7r™(a- 

^ m+1 

(to + 1)! 



da 



cxp(— Aa)da 



(50) 



(51) 



toG 



AR' 



(to+ 1)! 



(to+ 1)! 



E 



(m- 1)! 



(52) 



m + 1 



exp( — -)- 

TO + 1 

m — 1 



exp(-AAn) J2 ^ ( 



i?2 



m + 1 



It is worth noting that a relay set may contain the same loca- 
tion for different relays. This can be interpreted as employing 
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the same node in different frequency bands for forwarding, 
although this is not common in practical routing. We notice 
that these sets form (™) hyperplanes in the 2m dimensional 
hyperspace, which are of measure 0. Hence, even if we require 
distinct relays and take the integral over feasible regions, we 
will stiU get the same results. 



Appendix D 
Proof of Lemma[2] 

We proceed in a similar spirit as in the proof of Lemma 
[l] Define n = (rii, ■ ■ • ,n,n+i)'^ and k = (fci, • • • ,km+i f- 
When the ith hop is of distance and ki attempts are em- 
ployed in the ith hop, the probability for successful reception 
is given by 



g iri,ki, ■ • • , rm+i,kjn+i) 

m+l 



m+l 

n 



G"' exp {^XtKn^r^) 



n: l-<n-<k 



<k 

/ ?n+l \ m+] 

exp -XtK E n^r^ H 



i=l 



(53) 



where . Therefore, we redefine X^um, Ysum to 

be 



ni \ xi + 



R 



)2 m-l 
i=i 



R 



The Cauchy-Schwartz inequality indicates 



f m+l \ / m+l 



E 

, 1=1 



m+l 



i=l 



^ X* > ^ 

~^ ^sum — „„i+i I ■ 



(54) 



Therefore, the Lebesgue measure of the ellipsoid d'^^{Z„i)<a 
can be calculated as 



TT a — 



VM* 



i:t=y 1/". 



midet (A*J 



(55) 



where A*^ is the canonical form corresponding to y^*,^ and 
can be written as 

(ni+n2 -n2 ... \ 

— 712 n2 + —713 • ■ ■ 

—773 773 + 774 ... 

V ... 77,„ + 7l,„+l / 

By Laplace expansion of the determinant, we get 

det(A;„) = (77,„ + 77„.+i) dct(A;„_i) - 77^ det(A,*„_2)> 

/m+l \ /m+l ^ \ 

^ det(A:„) = ( n "0 ( E - 1 ' (56) 

which follows by induction. Define a = X^S^ "^i^f ■ Taking 
an integral over different isosurfaces yields 

dy„*(a) 



hn = 



da 



Gxp (— Aa) da 



-A 



AH 

l+l 



(57) 



By combining ( 53 1 and pT\ , we can derive the average 
number of relay sets when retransmitting ki times in the j* 
hop as 

Ek(A^m) 

=A"(l-7)'" E (-l)™'*'^- 

n: l^n^k 

m+l / s 



n: l-<n-<k 



m+l ^ , X 7r"'7)T,exp J -^^^sTl-^/j.^ 



n 



.=1 ^ w A'« (nr=r (Er=r 

(-l)'"+V™(l-7)^ 



^m^r7 



E 



-G^^j=i 



n: l-<n-<k 



i=l 



(58) 



In addition, we can impose a constraint on the maximum 
total number of attempts M without specifying the number 
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of transmissions for each hop. For a typical relay set with 
the i* hop of distance r^, we denote by pi the success 
probability of hop i in any time slot. Among these M time 
slots, successful reception occurs when there exists m+1 slots 
< i < m + 1) that satisfy: (1) transmission in the hop 
is successful at time tf, (2) for 1 < i < j < m + 1, we have 
1 < U < tj < M. We apply a greedy approach to search for 
all possible scenarios that allow successful reception, which 
can be determined by the smallest t = {ti,--- ,tm+i) that 
satisfies the above two requirements. By "smallest" we mean 
there is no t ^ t that meets the requirement. This is identical 
to finding the interval (ii — 1, ^2 — ^i — 1, • • ■ , ^m+i ~ t,n ~ 1), 
or equivalently, finding a vector j = (ji, • • • ,jm+i) such that 
j ^ and • 1 < M — m—1. Hence, the success probability 
with m + 1 hop routing can be calculated as 



immediately derive 
A 



1 



< 



m In '^^^^ ^'^ 



InG — ln(m + 1) — In In ■ 



(62) 

which yields ( p7| . Furthermore, in order to make the capacity 
well-defined, i.e. Tm{e) > 0, we will have the constraint for 
outage probability stated in the corollary. 

The derivation in a best-effort setting is exactly the same. 



' m+l } 

m+1 I m+l 

«=1 j^ l<J\/-m-l «=1 

/m+1 \ 

- Up] e Ei-D'^^-^n 



(59) 



m+1 / . \ m+1 



,4=1 / j^.l<Af-m-l0^k^j 



n 



1 ^ i=l 
m+1 / . \ m+1 

E i-^r- E nf;,)n^^"^ 



k-" •l<Af-m-l 
k)-0 



j^k 1=1 

j^l<Af-m-l 



1=1 



m+1 / . \ 771+1 

= E (-i)"^"^ E ^(^;:0^^^• (^o) 

k^-l<A/ j^ l<Af '=1 »=1 



Therefore, we can obtain 



m+1 



k" IKM 
k>-l 



j^-l<A/ ' = 1 



where Ek(-/Vm) is given in (58 i 



Appendix E 
Proof of Corollary[2] 

Let Dm — > oo, then the integral part in (j52]l becomes a 
Gamma function: 

lim / x™~ie^="da; = r(m) = (m - 1)! 



By setting the outage probability to be e, we can simplify ( 19 1 



- ^1 7™is:™(TO+l) m + 1 J 

Notice that the effective spatial density is Xj/{m + 1) and 
that e is monotonically increasing with respect to A, we can 
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